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Abstract. New observations and derived chemical abundances are reported for a sample of 57 bulge planetary 
nebulae (PN). Together with our previous results, a total of over a hundred objects have been analyzed, which 
constitute one of the largest samples of bulge nebulae studied under homogeneous conditions, including equipment 
and reduction procedures. In general, our data show a good agreement with some recent results in the literature, 
in the sense that the average abundances of bulge PN are similar to those from disk objects, however showing a 
larger dispersion. 
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1. Introduction 

■ In the past few years, many papers have been published 
dealing with the kinematics and abundances of the galac- 
tic bulge. Most of these works on bulge abundances are 

| concerned with heavy elements produced by supernovae, 
so that light elements such as helium and nitrogen have 
had a smaller share of attention. 

Planetary nebulae (PN) constitute an important tool 
in the study of the chemical evolution of the bulge, pro- 
viding accurate determinations of the abundances of light 
elements produced by progenitor stars of different masses. 
In fact, PN offer the possibility of studying both light ele- 
ments produced in low mass stars, such as He and N, and 
also heavier elements which result from the nucleosynthe- 
sis of large mass stars, such as oxygen, sulfur and neon, 
which are present in the interstellar medium at the stellar 
progenitor formation epoch. 

Previous determinations of chemical abundances of 
bulge PN have shown that, on average, these objects have 
abundances similar to disk planetary nebulae (Escudero & 
Costa 2001, Liu et al. 2001, Cuisinier et al. 2000, Costa & 
Maciel 1999, Ratag et al. 1997). However, the total num- 
ber of bulge PN with accurate abundances is still small; 
furthermore, a correlation between the chemical abun- 
dances and bulge kinematics is still to be determined, in 
contrast with the observed properties of the galactic disk. 
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In the present paper, we report new observations and 
derive chemical abundances for a sample of 57 bulge PN. 
These results are compared with previous data from our 
own group and other groups as well, both regarding the 
galactic bulge and other galactic systems, such as the 
galactic disk and halo. 

In section 2 we present our new observations for a sam- 
ple of bulge PN and comment on the reduction procedure. 
In section 3 we present our method to derive the chemical 
abundances. A discussion of the main uncertainties of the 
physical parameters is given in section 4. A comparison of 
our results with published data is given in section 5, and 
in section 6 we discuss our results an present our main 
conclusions. 

2. Observations and data reduction 

2.1. Observations 

In this work we present results for 57 new PNe, extend- 
ing and completing a previous paper (Escudero & Costa 
2001), in which we have presented our results for a sample 
of 45 objects toward the galactic bulge. 

The objects were selected according to their distances, 
taken from the literature (Schneider & Buckley 1996, van 
de Steene & Zijlstra 1995 and Zhang 1995), adopting as 
probable bulge objects those with heliocentric distances 
greater than 5 kpc. In spite of the sometimes larger un- 
certainties for these distances, this is still the best ho- 
mogeneous criterion to define the sample. Other criteria 
like diameters, H/3 or radio fluxes are even less reliable. 
It should be noted that in this paper we did not use the 
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Fig. 1. Distance distribution for our sample a (objects 
with heliocentric distances greater than 5 kpc) and objects 
from the literature. 

same criterion used by Escudero & Costa (2001), where 
all objects towards the bulge from Bcaulieu et al. (1999) 
and Kohoutek (1994) were observed. These objects do not 
have distance determinations, but as most of them have 
diameters smaller than 20" , they should be at or near the 
bulge according to the criteria of Gathier et al. (1983). 

Adopted distances are listed in table 3. We divided our 
sample in two samples: those with heliocentric distances 
greater than 5 kpc (Sample 1), and those with heliocen- 
tric distancies smaller than 5 kpc or unknown (Sample 2). 
Figure 1 displays the distance distribution for our Sample 
1 objects, combined with objects chosen from the litera- 
ture (see section 5 for detailed references). Adopting these 
distances and the criteria mentioned above we would ex- 
pect some 80 % of the sample to be in the galactic bulge, a 
fraction similar to the estimates given by Pottasch (1990), 
to which the reader is refered for a detailed discussion on 
the characterization of galactic bulge PN. 

The program objects were observed with two tele- 
scopes: 1.60m LNA (Laboratorio Nacional de Astroffsica, 
Brasopolis - Brasil), and 1.52m ESO (European Southern 
Observatory). Each object was observed at least twice in 
order to secure good final spectra. We selected the expo- 
sure time for each object in order to have a good S/N in 
the region of the H7 line, due to the importance of the 
[OIIIJA4363 line in the temperature diagnostics. 

For all objects, we used the long slit of 2 arcsecs width. 
At LNA we used a 300 1/mm grid with 4.4 A/pix disper- 
sion, and at ESO the used grid was 600 1/mm with 2.4 
A/pix dispersion. Each night at least three spectrophoto- 
metric standard stars were observed to improve the flux 
calibration. 

Table 1 list the PN G designation (first column), the 
usual names of the observed objects (column 2), equatorial 
coordinates for epoch 2000 (columns 3 and 4), date of 
observation (column 5), observatory (column 6) and the 
exposure time in seconds (last column). For some objects, 



spectra were taken with different exposure times: a short 
exposure to measure the Balmer series lines (and therefore 
to derive interstellar extinction), and a longer exposure 
(usually with Ha saturated) to measure the fluxes of lower 
intensity lines. 

Table 2 (available electronically) display the line fluxes 
of the reddened lines. For the cases where the line blend 
[OII]AA7319,7330 was impossible do separate, their inten- 
sities are added in the [OII]A7319 column. 

Figure 2 shows typical spectra from our sample. The 
upper panel shows M 3-26, an object with well defined 
lines, and the other shows KFL 4, a poorer object, ob- 
served with lower resolution. 



2.2. Interstellar Extinction 

The application of the extinction curve of Fitzpatrick 
(1999) had better results than the use of the curve of 
Cardelli et al. (1989) (CCM), when we compared the pre- 
dicted H7/H/3=0.47, (adopting Osterbrock's case B) with 
the calculated values after interstellar extinction correc- 
tion. A small difference between the two curves in the H7 
and Ha regions (sec fig. 6 of Fitzpatrick 1999) results in a 
significant discrepancy with respect to the recombination 
values when we use the CCM curve. Using the Fitzpatrick 
(1999) curve this difference is considerably reduced. These 
results are displayed in figure 3, where both distributions 
are plotted. Is easy to see that the distribution of H7/H/3 
ratios calculated with the Fitzpatrick curve is narrower 
and that its center is closer to the theoretical value of 
0.47. 

In order to bring the value of H7 closer to the recombi- 
nation value using the CCM curve, it would be necessary 
to use a value of the ratio of total to selective absorption 
Ry higher than 6.0, which is extremely large and inconsis- 
tent with the value determined by Stasihska et al. (1992). 
The main consequence for the poor interstellar correction 
in the Hy region is the poor determination of the elec- 
tron temperature for ++ . High 0[III]A4363 values would 
result in overestimated temperatures and low abundance 
values, as shown by other authors (Koppcn et al. 1991, 
Stasihska et al. 1998) who expected that some of the ob- 
jects lacking 0[III]A4363 might be the most metal rich. 

In this work we have used the Fitzpatrick (1999) ex- 
tinction curve, deriving E(B-V) from the observed Balmer 
ratio Ha/H/3 and adopting the theoretical value Ha/H/3 
= 2.85, with Ry = 3.1. 

Figure 4 displays the distribution of E(B-V) derived 
for Sample 1 (see section 2.1) compared with data from 
the literature and from Escudero & Costa (2001). It can be 
seen that the distribution for the present data agrees well 
with our previous work, and is in reasonable agreement 
with data from the literature, which are on average nearer 
the galactic center. 
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PN G Name RA(2000) DEC (2000) Date of Obs. Place Exp. Time (s) 



000.4-02.9 


M 3-19 


17 


58 


18 


- 30 


00 


10 


19-may-01 


ESO 


2x 


900 


000.4-01.9 


M 2-20 


17 


54 


25 


- 29 


36 


09 


19-may-01 


ESO 


lx 


600 


000.7+04.7 


H 2-11 


17 


29 


26 


- 25 


49 


06 


ll-jun-02 


LNA 


2x 


1800 


000.9-02.0 


Bl 3-13 


17 


56 


02 


- 29 


11 


15 


22-may/ll-jul-01 


ESO 


2x 


900, 2x 2400 


001.2+02.1 


He 2-262 


17 


40 


15 


- 26 


43 


21 


19-may-01 


ESO 


2x 


1200 


001.7+05.7 


H 1-14 


17 


28 


02 


- 24 


25 


28 


12-jul-01 


ESO 


2x 


1800 


002.0-02.0 


H 1-45 


17 


58 


23 


- 28 


14 


54 


19-may-01 


ESO 


3x 


600 


002.2-06.3 


H 1-63 


18 


16 


19 


- 30 


07 


36 


13-jul-01 


ESO 


2x 


1200 


002.2-02.5 


KFL 2 


18 


01 


00 


- 28 


16 


19 


22-may-01 


ESO 


2x 


1800 


002.6+08.2 


H 1-11 


17 


21 


18 


- 22 


18 


36 


14-jul-01 


ESO 


2x 


1200 


003.0-02.6 


KFL 4 


18 


02 


52 


- 27 


41 


02 


22-may/15-jul-01 


ESO 


2x 


1200, lx 3600, lx 1800 


003.3-07.5 


KFL 19 


18 


23 


09 


- 29 


43 


25 


12-jul-01 


ESO 


2x 


1800 


003.6-02.3 


M 2-26 


18 


03 


12 


- 26 


58 


31 


22-may-01 


ESO 


2x 


1200 


004.0-05.8 


Pe 1-12 


18 


17 


43 


- 28 


17 


24 


13-jul-01 


ESO 


lx 


1200, lx 1800 


004.8-05.0 


M 3-26 


18 


16 


09 


- 27 


14 


50 


12-jul-01 


ESO 


2x 


900 


004.8+02.0 


H 2-25 


17 


49 


02 


- 23 


42 


57 


19/21-may-01 


ESO 


lx 


900, lx 1200 


005.5+02.7 


H 1-34 


17 


48 


08 


- 22 


46 


48 


21-may-01 


ESO 


2x 


1200 


005.7-05.3 


M 2-38 


18 


19 


26 


- 26 


35 


20 


13-jul-01 


ESO 


2x 


120 


006.4+02.0 


M 1-31 


17 


52 


41 


- 22 


21 


57 


19-may-01 


ESO 


lx 


600, 2x 300 


006.8+02.3 


Th 4-7 


17 


52 


21 


- 21 


51 


12 


19-may-01 


ESO 


2x 


900 


007.0-06.8 


Vy 2-1 


18 


27 


59 


- 26 


06 


48 


13-jun-01 


LNA 


2x 


300 


008.1-04.7 


M 2-39 


18 


22 


02 


- 24 


10 


41 


ll-jun-02 


LNA 


2x 


300 


008.2+06.8 


He 2-260 


17 


38 


57 


- 18 


17 


35 


12-jul-01 


ESO 


lx 


1200, lx 600 


008.8+05.2 


Th 4-2 


17 


46 


09 


- 18 


39 


33 


12-jul-01 


ESO 


2x 


1800 


009.4-09.8 


M 3-32 


18 


44 


43 


- 25 


21 


34 


12-jul-01 


ESO 


2x 


900 


009.6-10.6 


M 3-33 


18 


48 


13 


- 25 


28 


56 


13-jul-01 


ESO 


2x 


600 


009.8-04.6 


H 1-67 


18 


25 


06 


- 22 


34 


51 


ll-jun-02 


LNA 


2x 


900 


010.7-06.7 


Pe 1-13 


18 


34 


53 


- 22 


43 


18 


13-jun-01 


LNA 


lx 


600, 2x 1200 


011.0+06.2 


M 2-15 


17 


46 


54 


- 16 


17 


27 


12-jul-01 


ESO 


2x 


900 


011.1+11.5 


M 2-13 


17 


28 


34 


- 13 


26 


17 


13-jul-01 


ESO 


2x 


1200 


011.3-09.4 


H 2-48 


18 


46 


35 


- 23 


26 


41 


10-jun-02 


LNA 


2x 


30, 2x 300 


013.4-03.9 


M 1-48 


18 


29 


30 


- 19 


06 


51 


10-jun-02 


LNA 


2x 


1200 


013.7-10.6 


Y-C 2-32 


18 


55 


31 


- 21 


49 


39 


13-jul-01 


ESO 


2x 


600 


018.6-02.2 


M 3-54 


18 


33 


04 


- 13 


44 


21 


ll-jun-02 


LNA 


2x 


1200 


018.9+03.6 


M 4-8 


18 


12 


06 


- 10 


43 


05 


ll-jun-02 


LNA 


2x 


1800 


018.9+04.1 


M 3-52 


18 


10 


30 


- 10 


28 


60 


ll-jun-02 


LNA 


2x 


1800 


352.0-04.6 


H 1-30 


17 


45 


08 


- 38 


08 


55 


10-jun-02 


LNA 


2x 


1200 


352.1+05.1 


M 2-8 


17 


05 


31 


- 32 


32 


08 


13-jun-01 


LNA 


2x 


800 


352.8-00.2 


H 1-13 


17 


28 


28 


- 35 


07 


32 


22-may-01 


ESO 


2x 


1200 


353.7+06.3 


M 2-7 


17 


05 


13 


- 30 


32 


14 


ll-jul-01 


ESO 


2x 


900 


353.8-01.2 


K 6-3 


17 


35 


27 


- 34 


47 


42 


14-jul-01 


ESO 


lx 


2400, lx 3600 


354.6-01.7 


K 6-5 


17 


38 


54 


- 34 


27 


39 


15-jul-01 


ESO 


2x 


3600 


356.1+02.7 


Th 3-13 


17 


25 


19 


- 30 


40 


44 


20-may-01 


ESO 


2x 


1200 


356.3-00.3 


Th 3-34 


17 


37 


45 


- 32 


15 


29 


14-jul-01 


ESO 


lx 


3600, lx 2400 


356.7-06.4 


H 1-51 


18 


04 


30 


- 34 


58 


00 


15-jul-01 


ESO 


2x 


1800 


356.7-04.8 


H 1-41 


17 


57 


19 


- 34 


09 


50 


10-jun-02 


LNA 


2x 


1200 


357.0+02.4 


M 4-4 


17 


28 


50 


- 30 


07 


46 


21-may-01 


LNA 


2x 


1200 


357.1-04.7 


H 1-43 


17 


58 


15 


- 33 


47 


39 


ll-jun-02 


LNA 


2x 


300, 2x 1200 


357.2-04.5 


H 1-42 


17 


57 


25 


- 33 


35 


44 


ll-jun-02 


LNA 


2x 


150 


ooo.o-uz.o 


Al 9 D 


1 7 


u i 


-to 




(Vi 




99 mav 01 

zz-iiiay-ur 




9 Y 
zx 


i son 


358.5-04.2 


H 1-46 


17 


59 


02 


- 32 


21 


44 


ll-jun-02 


LNA 


2x 


150, lx 600 


358.5-02.5 


M 4-7 


17 


51 


45 


- 31 


36 


00 


22-may-01 


ESO 


2x 


1800 


358.6+01.8 


M 4-6 


17 


35 


14 


- 29 


03 


11 


13-jul-01 


ESO 


2x 


1800 


358.8+04.0 


Th 3-15 


17 


27 


09 


- 27 


44 


20 


ll-jun-02 


LNA 


2x 


1800 


359.8+05.6 


M 2-12 


17 


24 


01 


- 25 


59 


23 


ll-jul-01 


ESO 


2x 


900 


359.8+06.9 


M 3-37 


17 


19 


13 


- 25 


17 


15 


12-jul-01 


ESO 


2x 


1800 


359.9+05.1 


M 3-9 


17 


25 


43 


- 26 


11 


54 


13-jul-01 


ESO 


2x 


1800 
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Table 2. Line fluxes (available electronically) 



Line 


Al 2-0 


Bl 3-13 


Hl-11 


Hl-13 


Hl-14 


Hl-30 


Hl-34 


Hl-41 


Hl-42 ... 


[OII]A3726+29 


179.3 




7.4 




48.0 


83.1 


57.7 


31.7 


18.5 ... 


[NeIII]A3869 


124.5 


49.3 


76.6 


74.8 


100.8 


112.9 




59.4 


88.5 ... 




o 0.3 



x 

Z3 



0.2 - 



0.1 



0.0 




6400 6500 6600 6700 6800 





6400 6500 6600 6700 6800 



_L 



1 1 

M 3-26 



KFL4 



_L 



3500 4000 4500 5000 5500 6000 6500 7000 7500 

Wavelength [A] 

Fig. 2. Typical spectra from our sample. The upper panel shows M 3-26, with well defined lines, and the lower panel 
shows KFL 4, a poorer object. 



3. Determination of chemical abundances 

3.1. Physical parameters 

The quality of the physical parameters is crucial for a good 
derivation of chemical abundances, as line emissivities of 
the ions are extremely sensitive to the electron tempera- 
ture and density of the nebulae. 

For the determination of electron temperatures we 
used the following line ratios for oxygen and nitrogen: 
[OIII]AA4363,5007 and [NII]AA5755,6584, from which wc 
can obtain estimates of the electron temperature for high 
and low ionization regions. The electron density was de- 
rived from the sulfur line ratio [SII]AA6716,6731, giving 
therefore an average value for the whole nebula. Excitation 



classes were also defined for the object following the 
Webster (1988) criterium. For each nebula the H/3 sur- 
face flux was calculated using the derived H/3 flux and the 
fraction of the nebula within the slit. 



In table 3 we give the observed (column 2) H/3 in- 
tensities and corrected by interstellar extinction (in erg 
cm~ 2 s _1 ) (column 3), diameters in arcsecs from Acker et 
al. (1991) (column 4), H/3 surface flux (in erg cm~ 2 s _1 
arcsec -2 ), the interstellar extinction E(B-V) (column 6), 
the electron density in 10 3 cm~ 3 (column 7), the electron 
temperatures in 10 4 K (columns 8 and 9), excitation class 
(column 10) and heliocentric distances in kpc (column 11). 
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by the Fitzpatrick (1999) (hatched) and Cardelli et al. 



This Work - Sample 1 
Escudero & Costa 2001 
Literature 




Fig. 4. Distribution of E(B-V) values for our sample. Also 
included are the distribution of Escudero & Costa (2001) 
and data from the literature. 



3.2. Ionic abundances 

In order to calculate abundances for the ions N + , O , 
0+, 0++, S+, S++, Ar++, Ne++, we have used the fits 
by Alexander & Balick (1997) that show a good agree- 
ment with ionic abundances obtained from statistical bal- 
ance equations. For the helium ions (Hc + , Hc ++ ) we 
used the recombination coefficients from Pequignot et al. 
(1991), with the coefficients for He + collisional excitation 
correction derived by Kingdon & Fcrland (1995). Only 
for Cl ++ , not present in the Alexander & Balick (1997) 
tables, we used the statistical balance equation for the 
[ClIII]AA5518+5538 lines. 

In particular, the + abundance can be derived from 
two pairs of lines: A3727+29 and A7319+30. There is how- 
ever a systematic difference between the resulting values, 
with some tendency to higher abundances for the red pair. 
These differences are shown in figure 2. In this figure, the 



dashed line is a linear fit for our sample. To explain this 
difference, observed also by other authors, we raised four 
alternatives: interstellar extinction, data reduction, elec- 
tron temperature and density variations, and ++ recom- 
bination. 

- Interstellar Extinction: one of the factors that could 
cause this difference is a poor interstellar extinction cor- 
rection. To check this, we corrected all spectra with two 
interstellar extinction laws (Fitzpatrick 1999 and CCM 
1989) using two extreme values for Ry: 2.1 e 5.9. The re- 
sults do not imply a significant variation of the ++ abun- 
dance in any case. The only possibility to invoke extinction 
to explain such a difference would be a new, totally dis- 
tinct interstellar extinction law, which seems unrealistic, 
since most of the extinction comes from the disk 

- Data reduction: as the two line pairs are in opposite 
regions of the spectrum, poor flux calibration could affect 
the final data. However, it is unlikely that this effect could 
affect all data. Furthermore, this effect appears also in 
other, independent works (Stasinska et al. 1998). 

- Temperature and density variations: this can be 
an important factor in abundance calculations (Viegas 
& Clegg 1994). Regions with high density in the inner 
part of a planetary nebula contribute to the intensity of 
A7319+30, but not to A3727 (Mathis et al. 1998). With 
respect to the density, as the intensity of the [SII] line 
relative to H/3 depends on the density (Alexader & Balick 
1997), this value will represent preferentially lower density 
regions, and this effect would be more evident for bulge 
objects due to their smaller angular sizes. For disk objects, 
the slit would cover only the central part of a planetary 
nebula and the effect is not as important as for bulge neb- 
ulae. 

- Recombination of ++ : For nebulae with low density 
and electron temperature, dielectronic recombination of 
++ may play an important role (Rubin 1986, Aller & 
Keyes 1987), leading to differences in the derived ionic 
abundances. 

Figure 5 displays a comparison between the + ionic 
abundance derived from the blue (3726+29) and the 
red (7320+30)lines, both for our sample and that from 
Stasinska et al. (1998). In the figure, the continuous line 
is a y=x plot and the dashed line is a linear fit for all 
the points. As can be seen, there is a small discrepacy 
between both determinations, with a tendency for higher 
abundances when the red lines are used. The same effect 
appears for both samples. As the final explanation is still 
an open issue, we adopted the blue pair to derive the + 
abundance like most of the other works in the literature: 



log(Ot r2r /R + ) = (-0.1±0.3)+(1.06±0.08)log(O£, 25 /H 4 



The derivation of sulfur abundance requires both S + 
and S ++ ions. For those objects where S ++ was not avail- 
able we adopted the same technique used by Kingsburgh 
& Barlow (1994), deriving the S ++ abundance from the 
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Fig. 5. A comparison of the + ionic abundance from the 
blue (3727+29) and red (7319+30) lines, respectively. 

relation between ratios S ++ /S + and ++ /0 + . This rela- 
tion was derived from our data, as can be seen in figure 6. 
In the figure, filled and open circles represent the points 
included in the fit, and crosses are abundances poorly de- 
termined, which were not included. The dotted line is a 
y=x relation and the continuous line is the best fit repre- 
sented by the relation: 

log{S ++ /S+) = (0.21 ± 0.08) + (0.73 ± 0.06)log(O ++ /O+) 



Since it is not possible to get the chemical abundances 
of all ions of a given element, we use the ICF method 
(Ionization Corrections Factors) to derive the elemental 
abundances, as in our previous works. For details, see 
Escudero & Costa (2001), Costa ct al. (1996) and refer- 
ences therein. 

The adopted ICF for oxygen and nitrogen are those 
suggested by Peimbert & Torres-Peimbcrt (1977); for sul- 
fur and neon we used the ICF provided by Kingsburgh & 
Barlow (1994), except for those objects for which S ++ is 
not available; in this case we adopted the value derived by 
the relation given in the last section. For argon we used 
the ICF provided by Freitas Pacheco et al. (1993). These 
ICFs are listed below. 



He 
H 



He+ He++ 
+ ■ 



H+ 



H+ 



9. -(911 0+ \ Hc 



N 
H 



N+ O 
H+ 0+ 
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Fig. 6. Relation between S++/S+ and 0++/0+ 



3.3. Elemental abundances 

The final abundances are listed in table 4. The objects 
separated by an horizontal line in the final of table 4 and 
also in table 3 are those from our Sample 2 (see section 
2.1). In the discussion on the abundance distribution (see 
section 5) we do not include these objects. 



4. Errors in the physical parameters and chemical 
abundances 

In Escudero & Costa (2001) typical errors were derived 
through a Monte Carlo simulation, based on the line in- 
tensities and their errors, which requires the same expo- 
sure time for all measurements. This is not the case for 
the present work, so we examined the effect of these un- 
certainties on the abundances. 

From Escudero & Costa (2001), we have estimated av- 
erage errors to the abundances of 0.023, 0.18, 0.17, 0.24 
and 0.19, respectively for He, N, O, S and Ar. However, it 
should be remembered that errors are highly affected by 
interstellar extinction, in the sense that objects that have 
higher extinction should have greater errors in the abun- 
dances, on the average, depending also on the exposition 
time and the size of the telescope. Therefore, abundances 
of bulge objects, that have high extinction, typically have 
higher errors than disk objects, so that we can estimate 
the uncertainties in the abundances of bulge nebulae to 
be roughly 0.2 dex, apart from He/H. 

4.1. Physical parameters 

One of the main error sources in the chemical abundances 
are the uncertainties in electron temperatures and densi- 
ties, as discussed previously. Low S/N in weaker diagnostic 
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Table 4. Chemical Abundances 



Name 


He 


e(N) 


e(O) 


e(S) 


e(Ar) 


e(Ne) 


Al 2-0 


0.160 


8.43 


8.66 


6.57 


6.78 


7.95 


Bl 3-13 


0.105 


7.68 


8.91 


7.14 


6.54 


8.06 


H 1-11 


0.106 


8.02 


8.67 


6.87 


6.28 


7.95 


H 1-14 


0.103 


7.62 


8.27 


6.39 


5.66 


7.38 


H 1-30 


0.147 


8.78 


8.76 


6.66 


6.86 


8.13 


H 1-34 


0.013 


8.14 


8.45 


7.04 


6.34 




H 1-41 


0.123 


7.95 


8.64 


6.81 


6.22 


7.85 


H 1-43 


0.011 


8.44 


8.84 


6.80 






H 1-45 


0.099 


7.39 


7.52 


5.69 


4.97 


6.84 


H 1-46 


0.096 


7.93 


8.35 


6.95 


6.16 


7.40 


H 1-63 


0.059 


7.48 


8.06 


6.86 


6.09 


6.97 


H 1-67 


0.125 


8.50 


8.75 


7.09 


6.64 


8.00 


H 2-11 


0.151 


8.16 


8.29 


6.56 


6.64 




H 2-25 


0.081 


7.56 


8.26 


6.74 


6.70 




H 2-48 


0.048 


7.12 


7.93 


6.41 


6.03 




Hp 2-260 


0.009 


7.14 


8.18 


6.29 


6.97 




He 2-262 


0.118 


7.82 


8.48 


6.66 


6.24 


7.55 


KFL 2 


0.113 




8.24 




5.99 


7.33 


KFL 4 


0.099 




8.34 




5.71 


7.28 


KFL 19 


0.080 


7.38 


8.65 


6.93 


6.06 


7.64 


M 1-48 


0.143 


8.65 


8.84 


7.09 


6.88 


8.33 


M 2-7 


0.137 


7.91 


8.74 


6.61 


6.80 


7.91 


M 2-8 


0.223 


8.50 


8.54 


6.04 


6.40 


7.81 


M 2-12 


0.010 


7.96 


8.56 


6.60 


7.86 




M 2-13 


0.117 


8.22 


8.72 


6.96 


6.48 


8.07 


M 2-15 


0.128 


8.30 


8.61 


6.95 


6.38 


8.00 


M 2-20 


0.081 


8.07 


8.78 


6.83 


6.43 




M 2-26 


0.161 


8.77 


8.49 


6.88 


6.65 


7.53 


M 2-38 


0.152 


8.57 


8.58 


7.28 


6.73 


7.78 


M 2-39 


0.107 


6.94 


7.48 


5.45 


5.64 


6.76 


M 3-19 


0.101 


7.27 


8.00 




6.29 


7.22 


M 3-26 


0.137 


8.43 


8.40 


6.99 


6.36 


7.92 


M 3-32 


0.132 


8.27 


8.56 


6.95 


6.35 


7.87 


M 3-33 


0.072 


7.72 


8.62 


6.57 


6.00 


7.84 


M 3-37 


0.092 


8.34 


8.86 


6.72 


6.74 


8.34 


M 3-52 


0.129 


8.58 


8.79 


7.15 


6.91 




M 3-54 


0.122 


7.33 


8.67 


6.52 


6.24 


7.97 


M 4-4 


0.154 


8.74 


8.92 


7.25 


6.78 




M 4-6 


0.111 


8.47 


8.71 


6.98 


6.30 


8.01 


M 4-7 


0.102 


7.94 


8.33 


6.43 


5.79 


7.59 


M 4-8 


0.015 


7.24 


8.03 


6.24 


6.66 




Pe 1-12 


0.114 


7.90 


8.63 


5.41 


6.36 


7.82 


Pe 1-13 


0.099 


- 


9.03 


- 


6.66 


8.29 


Th 3-13 


0.089 


7.33 


7.61 


6.10 


5.98 


- 


Th 3-15 


0.107 


7.62 


8.06 


6.15 


5.89 


- 


Th 4-2 


0.180 


8.72 


8.50 


6.95 


6.75 


7.79 


Th 4-7 


0.072 


7.70 


8.35 


6.47 


5.90 


7.49 


Vy 2-1 


0.145 


9.04 


8.94 


6.23 


6.94 


8.10 


H 1-13 


0.123 


8.84 


8.99 


7.47 


6.67 


8.20 


H 1-42 


0.106 


8.27 


8.69 


6.97 


6.15 


7.94 


H 1 ^1 


o Dos 


8 97 






7 79 


y.oo 


K 6-3 


0.128 


7.98 


8.33 


6.40 


6.49 




K 6-5 


0.159 


9.12 


8.57 


7.32 


6.88 




M 1-31 


0.161 


9.02 


8.98 


7.69 


6.89 


8.28 


M 3-9 


0.059 


7.69 


8.68 


5.19 


6.34 


8.05 


Th 3-34 


0.193 


9.15 


8.14 


7.08 


6.43 




Y-C 2-32 


0.088 


7.59 


8.59 


6.61 


6.19 


7.81 



lines, as well as internal variations in the physical param- 
eters certainly affect the final abundances. As discussed 
in section 3.2, electron densities derived from [SII] can be 
underestimated since the emissivity of this line depends 
on density. 

4.2. Helium 

The main factors that affect the helium abundance are 
the uncertainties in densities and electron temperatures. 
Depending on these parameters, the correction due to col- 
lisional excitation can reach 50%. Therefore, errors in the 
physical parameters would considerably affect the chem- 
ical abundances. The existence of neutral helium would 
also have an influence on the derived abundances. This 
was tested on the basis of the criterion proposed by Torres- 
Peimbert & Peimbert (1977), according to which objects 
with log + /0 > -0.4 seem to have a substantial contribu- 
tion of neutral helium. For this sample, 11 objects fit this 
criterion, indicating that their derived helium abundance 
should be taken cautiously. 

4.3. Abundances from collisionally excited lines 

In addition to observational errors and those resulting 
from the inhomogeneity of the nebulae, the use of ioniza- 
tion correction factors (ICF) to estimate elemental abun- 
dances for the elements with collisionally excited lines will 
include additional error sources. 

For nitrogen, + is used as ICF, and small variations in 
its abundance will affect the nitrogen abundance. The use 
of 0^ 727 or Oy 325 ionic abundances will imply differences 
of 0.0 to 0.7 dex in the nitrogen abundance. 

The same effect occurs for sulfur, in whose ICF + 
is also present. Two observational effects affect the sul- 
fur abundance: the dependence of S + emissivity with den- 
sity, and the weak intensity of the line [SIII]A6312, which 
will result in a poor determination of the S ++ abundance. 
Therefore, the present sulfur abundances should be taken 
cautiously. 

For oxygen, argon and neon, whose lines are normally 
intense in planetary nebulae, the main error sources are 
the uncertainties in the physical parameters. 

5. Data from the literature 

To compare our data, we collected objects from 5 recent 
works in the literature: Stasihska et al. (1998), Cuisinier 
et al. (2000), Cuisinier et al. (1996), Samland et al. (1992) 
and Koppen et al. (1991). We use only those objects with 
distances already determined (Schneider & Buckley 1996, 
van de Steene & Zijlstra 1995 and Zhang 1995), and with 
heliocentric distances greater than 5 kpc. 

In figure 7 we show the distribution of chemical abun- 
dances of He, N and O from the literature and our data. 

The results from the first paper of this series are also 
shown (solid irregular lines). In order to homogenize the 
results, we redcrived the elemental abundances for those 
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data, using the same ionic abundances but adopting 
improvements used here: sulfur ICF from Kingsburg] 
Barlow (1994); + abundance derived preferentially fi 
the blue doublet, and S ++ derived acording the procec 
described in section 3.2, when necessary. 

These results have been obtained using the same 
struments and procedures as in the present paper, so t 
both samples are expected to present a large degrei 
homogeneity. This fact can be observed in fig. 7, as 
abundance distributions of both samples are similar. Sr 
differences observed at the low end of the He and N ab 
dances can be attributed to the incompleteness of the s; 
pies. It can be seen that our abundances show a good s 
ilarity to all other data. The nitrogen abundances, t 
reflect the progenitor mass, show a little more scattei 
from our objects than the data from literature. This 
expected to occur because our objects are more wic 
distributed in angular position than the literature obji 
(see Escudero & Costa 2001). 

6. Discussion and conclusions 

Our results are in agreement with those from other wc 
(see for example Escudero & Costa 2001, Cuisinier et 
2000), indicating that bulge PNe have mean abundai 
similar to those from disk objects, however displaj 
higher dispersion. 

Figure 8 displays the relation between excital 
classes and helium abundances (panel a) and oxygen ab 
dances (panel b). The data distribution in panel (a) sh 
no correlation between the helium abundances and 
excitation class; the only exception is that some obji 
with very low abundances are preferably associated v 
the lowest excitation classes, which possibly indicates 
presence of neutral helium for these objects. On the oi 
hand, the absence of correlation in panel (b) also indie; 
that our abundance analysis is unbiased with respeci 
the excitation of the nebulae. 

Figure 9 displays several abundance correlations 
our sample which are distance independent, combi 
with data from Escudero & Costa (2001) and the li 
ature. 

Oxygen and nitrogen are two key elements for the d 
nostics of PNe. Due to strong lines, oxygen abundances 
normally very well determined, giving the abundano 
the interstellar medium at the progenitor formation epc 
It was already noted that more massive, younger obji 
like type-I PNe progenitors have oxygen abundance lo 
than the Sun or type-IIa PNe progenitors (see Costa e1 
1996). The same effect appears in stellar evolution yi< 
(van den Hock & Groenewegen 1997, Marigo 2001), h 
ever for the majority of the lower mass stars, oxygen 
mains unchanged and can be used to trace the chem 
evolution of the interstellar medium. 

Nitrogen, on the other hand, has two possible orig 
it reflects the abundance of the interstellar medium at 
progenitor formation epoch, combined with the mate 
produced by nucleosynthesis during the evolution of 
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Fig. 7. Distribution of helium, nitrogen and oxygen abun- 
dances for this work (Sample 1), from Escudero & Costa 
(2001) and from selected sources in the recent literature 
(see text). 



PN progenitor star. Typical yields as a function of the 
stellar mass are given by van den Hoek & Groenewegen 
(1997) and Marigo (2001). 

Examining the O vs. N diagram displayed in fig. 9a, 
a positive correlation between O/H and N/H can be ob- 
served. As for the correlation between log(N/0) and O/H 
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(figure 9b), this is less evident. We can see that objects 
with high N/O can have higher or lower oxygen abun- 
dances. Those with higher oxygen abundances are similar 
to the disk type-I PNe, while those with lower abundances 
are similar to halo objects. On the other hand, objects 
with low N/O are usually oxygen-rich, which is consis- 
tent with the correlation shown in fig. 9a. Some of these 
objects may be type-Ill PNe in the inner disk (see also 
Maciel 1999) and with the abundances of metal rich bulge 
giants (Rich 1988). This dispersion in the abundances sug- 
gests that the chemical evolution of the bulge cannot be 
explained in a single formation scenario, so that a com- 
posite scenario is required instead. 

Another important diagram for the diagnostics of PNe 
abundances is He x log(N/0), which reflects the chemical 
enrichment during the stellar evolution of the progenitor. 
From figure 9c it can be seen that there is a reasonably 
well defined correlation between He/H and log(N/0), in 
the sense that helium-rich nebulae are are usually also ni- 
trogen rich, which probably reflects an increasing progen- 
itor mass. A small group of objects apparently have very 
low He abundances, while their N abundances are normal. 
This can be explained either by the fact that their progen- 
itor star had smaller masses, which characterizes them as 



type Ha nebulae, or their He abundances may have been 
underestimated owing to the presence of neutral helium. 

Figures 9d, 9e and 9f display the correlation between 
the alpha-elements; as they are produced mainly in type 
II supernovae their correlation reflects the chemical evo- 
lution of the interstellar medium from which the progeni- 
tor stars are formed. From these figures one can see that 
lower oxygen abundances imply also lower sulfur, argon 
and neon abundances. 
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